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Abstract 

A low-speecl wind-tunnel study of the flow about a 76/40-deg double-delta wing is described for angles of attack 
ranging from -10 to 25 deg and Reynolds numbers ranging from 0.5 to 1.5 Million. The study was conducted 
to provide data for the purpose of understanding the vortical flow behavior and for validating Computational 
Fluid Dynamics methods. Flow visualization tests have provided insight into the effect of the angle of attack and 
Reynolds number on the vortex-dominated flow both on and off of the surface of the double-delta wing. Fpper 
surface pressure recordings from pressure orifices and Pressure Sensitive Paint have provided data on the pressures 
induced by the vortices. Flowfield surveys were carried out at an angle of attack of 10 deg by using a thin 5-hole 
probe. Numerical solutions of the compressible thin-layer Navier-Stokes equations were conducted and compared 
to the experimental data. 


*The present study has been conducted as a cooperative research agreement between NASA Langley Research Center and the Naval 
Air Warfare Center. The first author was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NASI- 19480 while in residence at the Institute for Computer Applications in Science and Engineering (ICASE). NASA Langley 
Research Center, Hampton. VA 23681-0001. 
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1 Introduction 


St ,rake/ wing configurations are employed for several modern peace keeping aircraft and for design concepts of 
high-speed civil transports. The highly swept strake is an additional lifting surface above which a stable vortex is 
formed. This strake vortex persists over the wing and induces a strong outflow, which energizes the boundary layer 
on the upper surface of the wing and keeps it attached to a higher angle of attack than in the case of a wing w ithout 
a strake. Therefore, strakes increase the lift to drag efficiency of delta wings [1]. At the strake/ wing leading-edge 
kink or juncture, a w T ing vortex originates that increases in size and strength with distance downstream from the 
kink. The strake and wing vortices have the same direction of rotation and at high angles of attack their trajectories 
entwine due to interaction. A limit to these effects is reached once the angle of attack attains values where vortex 
breakdown occurs over the wing. 

Published experimental data on the structure of the subsonic flow over strake/wing configurations is available 
from several studies carried out in water and wind tunnels on models that are mostly double-delta wings. The 
results of these studies are used to validate numerical solutions. A summary of the experimental results and 
available numerical solutions is given by Rom in Ref. 2, and Verhaagen and Maseland in Ref. 3. The experiments 
indicate that the interaction process and breakdown of the strake and wing vortices depend on the angle of attack, 
the shape of the leading-edge cross section and the leading-edge kink angle. The latter angle is the difference 
between the leading-edge sweep of the strake and the leading-edge sweep of the wing. In general, the cross-over- 
point of the interacting strake and wing vortices moves downstream with decreasing angle of attack and increasing 
leading-edge kink angle and Reynolds number. At incidences and kink angles less than 10 deg, no interaction 
of the vortices occurs. The majority of the experimental and numerical studies were carried out on double-delta 
wings with strake leading-edge sweep angles between 76 and 80 deg and wing sweep angles near 60 deg. These 
sweep-angle combinations are required for high-speed aircraft configurations to achieve sufficient supersonic-cruise 
performance. 

For highly-maneuverable aircraft, strake/wing configurations are employed which feature smaller wing leading- 
edge sweep angles to increase wing area and low* speed performance. A 76/40-deg planform is representative of 
such a type of aircraft [4], Experimental data for double-delta wings that have a 40 deg swept wing are rather 
limited. Experimental studies by Graves, et al [1]. Brennenstuhl [5], Hebbar, et al [6], Thompson [7] and Verhaagen 
[8] provide qualitative insight into the effect of the angle of attack on the crossflow structure, the interaction and 
breakdown of the strake and wing vortices, and the forces and moments. However, no surface pressures or flowfield 
data are presented in these references. As far as numerical solutions are concerned, solutions of the Navier-Stokes 
equations for the flow about a double-delta wing with a 40 deg swept wing were computed by Kern [4], Hsu and 
Liu [9], and Ekaterinaris, et al [ 1 0] . The solutions predict the interaction and breakdown of the strake and wing 
vortices as well as the surface pressure levels reasonably well. For a more detailed validation of numerical solutions, 
additional experimental data are needed on the effect of angle of attack, sideslip and Reynolds number on the 
vortex flowfield and induced wing surface pressures. 

The objective of the present study was to extend the database for 76/40-deg double-delta wings and to compare 
the experimental data with numerical solutions. This report summarizes the results of the flow visualization tests 
and surface pressure measurements conducted at angles of attack ranging from -10 to 25 deg, and flowfield surveys 
conducted at an angle of attack of 10 deg. The experimental data for the latter angle of attack is compared with a 
recent numerical solution computed at the Naval Air Warfare Center (NAWC). Results of the present study have 
been published earlier in Ref. 11. 


2 Experimental Setup 

Wind Tunnel and Test Conditions 

The tests were carried out in the Basic Aerodynamics Research Tunnel (BART) of the NASA Langley Research 
Center (LARC). This is a subsonic, open-return wind tunnel with a test section 0.71 m high. 1.02 m wide and 
3.05 m long (Fig. 1). The test section consists of two 1.5-m cells, with large windows which facilitate the use 
of non-intrusive diagnostic and detection techniques. The airflow entering the test section is conditioned by a 
honeycomb, four anti-turbulence screens and an 11:1 contraction ratio. The variation in longitudinal turbulence 
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intensity ranges from approximately 0.03% at an air speed of 15 m/sec to 0.09% at an air speed of 48 m/sec 
Further details of the wind tunnel facility can be found in Ref. 12 by Sellers and Kjelgaard. 


Model 

The model was a stainless steel 76/40-deg swept double-delta wing designed and provided by NAWC (Fig. 2). 
The model had a chord length of 0.406 m, a span of 0.415 m and thickness of 9.53 mm. The leading-edge kink 
was situated at 52%. chord and the wing tip was cropped from the 84%> chord wise station. The sharp leading and 
trailing edges of the model were broken and rounded to a diameter of 0.13 mm (0.005 in.). The under surface was 
beveled at 20 deg over a width of 26.16 mm. The kink section was capable of allowing fillets of different leading 
edge geometry to be installed [4]. In the present study, only the baseline configuration with straight strake and 
wing leading edges was investigated. The left upper surface of the model was equipped with 48 pressure orifices 
of 0.25 mm (0.010 in.) diameter. Ten orifices were drilled in a spanwise row at 25% chord, while 19 orifices were 
drilled both at 75%) and 90%> chord. The pressure tubes exited the model at its lower surface. The model was 
supported in the tunnel by an internal-balance sting attached to its lower surface (Fig. 3). A fairing was fitted 
to the lower surface of the model to cover the balance block and the pressure tubes. The angle of attack of the 
model was controlled remotely and set with an accelerometer mounted on the sting at a short distance behind the 
trailing edge of the model. The maximum uncertainty in the angle of attack reading was 0.03 deg. At large angles 
of attack, the interaction between the flow over the model and the tunnel walls resulted in a lift interference effect 
that increased the angle of attack. The pressure field upstream of the sting model support induced an additional 
upflow. Based on data published in Ref. 13. it is estimated that these interference effects caused the effective angle 
of attack to be on the order of 2 deg larger than the geometric angle of attack. 

Vortex Flow Visualization 

A smoke/laserlight-sheet, technique was used to visualize the structure of the vortex flow and to study probe/vortex 
interference effects. A single column of vaporized propylene glycol was introduced upstream of the tunnel inlet and 
positioned such that it could be entrained into the vortices generated by the model. A 6- Watt Argon ion laser and 
a twin-mirrored galvanometer laserlight-sheet generator were used to illuminate the propylene glycol particles and 
thus the vortices. The vortices were viewed and recorded by two video cameras through a window on the top of 
the test section. The cameras were set at different viewing angles to capture the flow in planes both normal and 
parallel to the surface of the model. The laserlight-sheet technique was used at angles of attack ranging from 5 to 
22.5 deg and a Reynolds number (i?e) of mostly 0.5 million (M). 

During the flow visualization tests, the upper surface of the model was covered with black vinyl contact paper 
to prevent the propylene glycol from settling in the pressure orifices and to enhance the contrast. The covering 
was also used to overlay a grid on the upper surface of the model to mark the centerline and the stations from 10 
to 90%) chord, at intervals of 10%; chord. 

Surface Flow Visualization 

Surface flow visualization was used to determine the location of stagnation, attachment and separation lines 
in the boundary layer pattern on the upper surface of the model. This pattern was visualized using a mixture of 
titanium dioxide (Ti02) and kerosene. The mixture was applied to the model and the air speed in the test section 
was brought up to the test conditions as fast as possible. The development of the surface flow pattern was recorded 
on video tape. After the mixture on the model surface had dried completely, the contact paper was removed from 
the model and photographed. The surface flow pattern was investigated at angles of attack ranging from -10 to 
22.5 deg and Re = 0.5, 1.0 and 1.5 M . 

Surface Pressure Measurements 

The distribution of the upper-surface static-pressure coefficient Cp at the 25, 75 and 90%; chordwise stations 
was measured at angles of attack ranging from -5.0 to 25.0 deg and Re — 1.0 M. To study the effect of the Reynolds 
number on Cp, at 10 and 22.5 deg angle of attack (AO A) pressures were acquired at Re ranging from 0.5 to 1.5 M. 
While preparing the model for the surface pressure measurements, 2 of the 10 pressure ports at 25% chord were 
found to be unreliable. Therefore, these ports were not used during the testings. The remaining 46 pressure ports 
in the model were connected to transducers with a maximum range of +/-1.0 p.s.i. and an uncertainty of 0.2%). 
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The uncertainty in Cp is a function of the surface pressure and the dynamic pressure. For Re = 0.5, 1.0 and 1.5 M. 
the maximum absolute uncertainty in the measured Cp is estimated to be +/- 0.275, 0.079 and 0.031, respectively. 
The data remain uncorrected for tunnel wall and support interference effects. 

The Cp distribution on the entire upper surface of the model was visualized by applying Pressure Sensitive 
Paint (PSP) at AO A = 10 and 22.5 deg and Re = 1.0 and 1.5 M. The technique uses a surface coating containing 
probe molecules that luminesce when excited by an appropriate light source. Oxygen molecules interfere with this 
process and decrease the amount of luminescence. As a result, the luminescence of the paint varies as a function 
of the partial pressure of the oxygen. Therefore, the intensity of the luminescence can be related to the static 
pressure of the air at the coated surface. The paint was excited with eleven tungsten-halogen lamps outfitted with 
blue stacked filters. The photoluminescence of the paint was detected with a scientific grade digital camera with an 
amber filter. As part of the calibration of the paint, the ratio of images taken at two pressure conditions was ob- 
tained. This corresponds to a ” wind-on" condition where the pressures were unknown, and a v wind-ofF* condition 
where the model pressure distribution was constant and equal to atmospheric pressure. The offset resulting from 
the charge-coupled device (CCD) detector dark current and the noise levels of the camera was subtracted from the 
two images before the ratio was taken. The final paint calibration coefficients were then determined using a linear 
curve fit with three pressure ports on the model. Further details on the technique can be found in Refs. 14 and 15. 

Flowfield Surveys 

Quantitative data of the vortical flowfield over the wing was obtained by surveying the flow using a hemispher- 
ical, five-hole probe (5HP) having a diameter of 3.2 mm . The probe head had a pressure port on the centerline, 
while four ports, at angle approximately 40 deg off the probe centerline, were located symmetrically around the 
central pressure port. An iterative calibration procedure using modified equations for the flow around a sphere 
was used. This procedure is described in Ref. 16 by Kjelgaard. The probe data were reduced to the three velocity 
components and the static- and total-pressure coefficients. Using this calibration, the 5HP could be used to measure 
total flow angles with respect to the probe axis of up to 50 deg and flow velocities of up to 54 m/sec (180 ft/sec). 
The estimated measurement uncertainty for the total flow angle was 1 deg for flow angles less than 40 deg, and 2 
deg for flow angles ranging from 40 to 50 deg. The surveys were carried out at AO A — 10 deg in planes normal 
to the model upper surface at the 50, 62.5 and 75% chordwise stations. In order to keep the measured velocities 
within the calibration of the 5HP. flowfield data were acquired at a freestream velocity of 18 m/sec ( Re = 0.5 
M ). Earlier experiments indicated that the velocities in the vortex core can reach values of up to three times 
[3]. The 5HP was mounted on a five degree-of-freedom traverse (xyz- translation, pitch and yaw) with four degrees 
under computer control. The traverse enabled movement of the probe with a translation resolution of 10 Jim in 
the x . y and 2 directions, and a yaw resolution of 0.01 deg. For all of the survey planes, the probe was yawed to 
approximately align it with the vortex axis, thus reducing the effective total probe flow angle. 


3 Experimental Test Results 

Vortex Flow Visualization 

Fig. 4 shows pictures of the strake and wing vortices illuminated by a laserlight sheet at AO A = 10 deg. The 
pictures are photoprints of frames from the recorded video tape taken by the video camera that viewed the model 
from the front at an oblique angle. In the upper picture, the vortex crossflow at 75% chord was illuminated by 
the laserlight sheet orientated normal to the model. In the bottom picture, the laserlight sheet was orientated 
such that it passed through the vortex axes. The dark void along the axes of the vortices indicates a lack of 
light-scattering smoke particles in this region. The void could possibly be due to the high circumferential velocities 
near the edge of the subcore through which the smoke particles are centrifuged out of the core. In addition, the 
high axial velocities in the subcore tend to reduce the density of the smoke. From the picture it is evident that the 
strake and wing vortices hardly interact at AO A - 10 deg. Fig. 5 gives a planview of the vortex axes trajectories 
as well as the core bursting locations observed at the various angles of attack. The trajectories shown, suggest that 
the interaction between the strake and wing vortices increases with the angle of attack. This is due to the vortices 
increasing in size and strength. It can further be noted that the vortices burst before they cross over. At .40,4 
— 15 and 20 deg, the strake vortex burst when it is passing underneath the burst wing vortex. This may indicate 
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that the breakdown of the strake vortex is triggered by the burst of the wing vortex. The wing vortex cores hence 
have been observed to burst earlier than the strake vortex cores. This tendency has also been found in other wind 
and water tunnel experiments conducted with 76/40-deg swept and 80/40-deg swept double-delta wings [5,7 and 
8]. An exception is the water tunnel study of Hebbar [6], who observed the wing vortices to burst downstream of 
the breakdown point of the strake vortices. In Fig. 6, the range of observed breakdown locations is given for each 
of the angles of attack investigated in the present test. The breakdown point coordinate is denoted xbr and has 
been non-dimensionalized by the root chord length c. The breakdown location fluctuated over a distance of about 
5% root chord and moved upstream with increasing angle of attack. The location where the core started to expand 
was taken as the bursting point. Downstream of this location the core was observed to diffuse rapidly. Payne, et 
al, [17] distinguish a spiral- or bubble- type of vortex core breakdown. Unfortunately, the temporal resolution of the 
video recordings was not detailed enough to identify the type of breakdown occuring in the strake and vnng vortex 
cores. The effect of Reynolds number on the burst location of the strake vortex was investigated at AO A = 20 deg. 
The burst location was observed to be unaffected by a change in Re from 0.25 to 2.0 M. This confirms suggestions 
made by Erickson [18] that the onset of bursting is independent of Re. 

Fig. 7 shows a close-up of the strake vortex near the leading-edge kink at AO A = 20 deg. In the free shear layer 
small stationary vortices can be seen which emanate from the leading edge. The vortices are co-rotating and have 
the same direction of rotation as the primary vortex. This type of vortical flow structure in a free shear layer w^as 
for the first time noted by Squire[19]. His illustration has been reproduced here in Fig. 8. Detailed investigations 
of the vortical structure by Washburn and Visser[13] indicated that these vortices follow* streamlines in the flow as 
they follow a helical path around the primary vortex core. 

Surface Flow Visualization 

Fig. 9 shows pictures of the boundary layer flow pattern at AO A = -10, 0, 10 and 22.5 deg, as well as sketches 
indicating the characteristic features of the surface flow pattern at these angles. The boundary layer flow at negative 
and zero angle of attack was investigated at Re = 1.5 M only. At lower Re , the surface shear flow became too 
weak to generate a clear surface pattern. For positive angles of attack, the patterns were investigated at Re — 0.5. 
1.0 and 1.5 M. Only the patterns for Re — 1.0 M are shown here. The effect of Re on the flow T pattern will be 
discussed at the end of this section. 

At negative angles of attack, a stagnation line is discernable on the strake upper surface just inboard of the 
leading edges. Outboard of this line, the boundary layer is directed towards the leading edge where it separates 
and mixes with the boundary layer from the lower surface. Inboard of the stagnation line, the boundary layer 
remains attached and flows in the chord wise direction. It can be expected that there will also be an attachment 
line inboard of the wing leading edges. The resolution of the flow pattern in this region, however, w^as not detailed 
enough to discern such a line. Near the cropped leading edge, the surface flow' bends towards the tip. This flow 
mixes with the boundary layer flowing off the lower surface to form a free shear layer that rolls up into a tip vortex. 

At zero deg angle of attack, the boundary layer flow pattern shows outflow markings induced by weak strake 
and wing vortices. Although the upper surface is parallel to the freestream at this angle, support interference 
effects induce a small upflow angle that causes the flow to separate at the strake and wing leading edges. 

At positive angles of attack, the boundary layer flow pattern on the strake changes little with the angle of attack. 
The pattern is characterized by outflow markings induced by the strake vortices and clear secondary separation 
lines. The attachment lines of the primary vortex flow are outboard of the symmetry plane at low* angles of attack 
and reaches the wing centerline at AO A % 15 deg. The secondary separation lines are straight, giving evidence 
of an approximately conical strake vortex flow. The boundary layer underneath the strake vortex separates at a 
spanwdse position y/s « 0.70, indicating that the boundary layer is laminar ( s denotes local wing semispan ). In 
the case of a turbulent boundary layer, separation would occur more outboard at y/s % 0.85 [20]. The topology 
of the crossflow over the strake is shown in the left sketch of Fig. 10. The free shear layer of the secondary vortex 
attaches just inboard of the leading edge and induces an inward directed flow on the surface (see insert sketch 
AO A — 10 deg, Fig. 9). This flow* separates again at the tertiary separation line. Between the latter line and the 
secondary separation line a region of accumulated surface oil was visible in the experiment. At low Re, the surface 
shear in this region was too weak to generate a clear pattern. At high angles of attack and Re, the pattern showed 
markings of an outflow that is assumed to be induced by the tertiary vortex. The tertiary separation line w r as 
visible from near the apex to the leading edge kink. Downstream of the leading edge kink, this line disappears as 
the tertiary vortex is supposed to be entrained into the stronger wing vortex. 
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Downstream of the leading edge kink, the effects of the angle of attack on the surface pattern are much stronger. 
Up to .40.4 = 10 deg, the shear flow pattern induced by the shake and wing vortices and accompanying secondary 
separation lines remain individually distinguishable up to the trailing edge. Hie conjectured topology of the 
crossflow at .40.4 = 10 deg is shown in the right sketch of Fig. 10. It is characterized by a double-branched wing 
vortex core fed by vorticity from free shear layers that are connected to the wing leading edge and the strake vortex 
core. 

Beyond .40,4 = 10 deg, the strake vortex secondary separation lines vanish starting at the trailing edge, 
indicating that the secondary separation is suppressed. Both the strake and wing vortex cores break down over 
the wing. Just upstream of the bursting point of the strake vortex, the surface pattern gives evidence of the strong 
outward bending of this vortex. Downstream of the bursting point a wide region of t urbulent surface flow is visible. 
Outboard of the wing vortex secondary separation line elongated regions of accumulated and circulating surface 
fluid were observed. These regions tend to move towards the kink with increasing angle of attack. Near the wing 
tip a rather complex surface flow pattern is visible. This pattern is assumed to be generated by the bluff-body type 
wake flow' downstream of the burst strake vortex and by the wdng tip vortex. Near the tip a focus was observed in 
the surface flow pattern. This focus induced an upstream flow locally at the trailing edge. 

The effect of the Reynolds number on the surface pattern is illustrated in Fig. 11 for AO A = 10 and 22.5 deg. 
On the strake, the pattern changes little when Be increases from 0.5 to 1.5 M. There were no signs of boundary 
layer transition from laminar to turbulent. On the wing, the strake vortex secondary separation lines are clearly 
visible at ,40.4 = 10 deg when Be = 0.5 and 1.0 M. From about the 65 % chordwise position, the secondary sepa- 
ration lines at Be = 1.0 M are located outboard the lines at. Be = 0.5 M. This is a consequence of the transition 
of the boundary layer underneath the strake vortex from laminar to turbulent. The secondary separation lines 
have vanished at Be = 1.5 M, indicating that no secondary vortices are induced by the strake vortices at this Be. 
Another possibility may be that the secondary vortices have moved outboard and merged with the wing secondary 
vortices. At ,40.4 = 22.5 deg. on the wing no strake vortex secondary separation lines were visible at the Be 
investigated. Reynolds number effects were further noted on the spanwise location of the wing vortex secondary 
separation lines. These lines tend to bend outboard with increasing Be due to transition of the boundary layer 
underneath the wing vortex from laminar to turbulent. 


Surface Pressure Measurements 

Fig. 12 shows the effect of the angle of attack on the upper surface Cp at 25% chord for Re = 1.0M. The re- 
distribution has been plotted versus the wing lateral coordinate y, non-dimensionalized by the local wing semispan 
s. At AO A — 0 deg, a small suction peak is induced by the weak strake vortex near the leading edge. Inboard of 
the suction peak, flow incident to the wing causes a region of positive Cp. The angle of attack has little effect on 
the crossflow structure sketched in Fig. 10. The size and circulation of the strake vortex, however, increases with 
the angle of attack, resulting in a growing suction peak. Outboard of this peak the pressure curves rise slightly due 
to the influence of the secondary and tertiary vortices. The level of suction induced by these vortices is about the 
same as that induced by the primary vortex, confirming that the boundary layer underneath this vortex is laminar. 
In the case of a turbulent boundary layer, the suction induced by the primary vortex would be larger, while the 
opposite would be the case for the suction induced by the secondary vortex. There would probably be no tertiary 
vortex [21]. 

The effect of the angle of attack on the Cp-distribution at 75% chord, shown in Fig. 15. can be explained using 
Fig. 14, which illustrates the effect of the angle of attack on the topology of the crossflow at this station. The 
latter figure is based on the results of the flow visualization tests. 

At .40.4 — 0 deg, a weak strake vortex suction peak and stronger wing vortex suction peak can be detected in 
the Cp-distribution. As the angle of attack is increased to 5 deg, a slight increase can be noted in the strake vortex 
suction peak and a strong increase of the suction induced by the wing vortices. Outboard of both t he strake and the 
wing vortex suction peak, a suction peak is induced by the accompanying secondary vortex. Up to .40.4 = 5 deg, 
the strake and wdng vortices are small and not connected to each other. In the gap between these vortices flow is 
incident to the wing, thus causing a small region of positive Cp at about y/s ^ 0.6. This region of positive Cp has 
also been noted in the Cp data of an identical double-delta wing investigated at the Delft University of Technology 
(TUD) [22]. When the angle of attack is increased, the strake and wing vortex cores become interconnected. Up 
to AO A - 12.5 deg, the suction induced by the strake and wing primary and secondary vortices can be seen to 
increase. The increase in suction of the wing vortex secondary vortex is only small. This may be due to the 



transition of the boundary layer underneath the wing vortex from laminar to turbulent which reduces the influence 
of the secondary vortex. Beyond ,40.4 = 10 deg, the strake and wing vortices start to coil around each other. As a 
consequence, the strake vortex moves closer to the wing surface and more outboard, while the wing vortex moves 
upwards and inboard. The effect of these displacements is noticable in the Op-curves; up to AO A — 20 deg the 
suction peak of the strake vortex moves outboard and increases in size, while the suction peak of the wing vortex 
moves inboard and reduces rapidly. At .40.4 = 20 deg, the wing vortex is overhead the strake vortex and only one 
suction peak is induced on each wing half. At AO A = 25 deg, the pressure distribution is affected by the burst of 
the strake vortex core. As a result, the suction peak not only collapses, but also spreads in the span wise direction. 
The latter is caused by the radial expansion of the burst vortex core, through which the circulation previously 
contained in a narrow core is redistributed over a wider region. 

The effect of the angle of attack on the Cp measured at 90% chord is shown in Fig. 15. At this station, suction 
is less than at 75% chord for all angles of attack considered. Separate strake and wing vortex suction peaks were 
measured up to AO A — 12.5 deg. At 15 deg, these peaks have merged into a single peak, while at higher angles of 
attack the surface pressures become affected by vortex bursting. 

The effect of Reynolds number on Cp was investigated at AO A = 10 and 22.5 deg by acquiring data at Re 
ranging from 0.5 to 1.5 M. At both AO A, the effect of Re on the Cp distribution at 25% "chord is of the order of 
the measurement accuracy. At 75% chord, the upper plot of Fig. 16 show's that at AO A = 10 deg there is a clear 
effect in the region dominated by the strake vortex. Due to the transition of the boundary layer underneath the 
strake vortex from laminar to turbulent, the secondary vortex and its suction have become nil at Re — 1.5 M. This 
is in accordance with the surface flow pattern where the secondary separation lines have vanished at. this Re. As a 
result, the primary vortex is located closer to the wing surface, inducing a higher suction than at Re = 0.5 M. On 
the wing vortex side, transition of the boundary layer results in a smaller and more outboard located secondary- 
vortex. As a consequence, the primary wing vortex is expected to move slightly outboard as well, explaining the 
slight outboard shift of the wing vortex suction peak that can be noted when Re is increased from 0.5 to 1.5M. 
The lower plot of Fig. 16 shows that at AO A = 22.5 deg Re has a similar effect on the suction induced by the 
strake primary and secondary vortices as noted at AO A = 10 deg. 


4 Computations 

About, the double-delta wing a spherical C-0 topology grid was created using a three-dimensional hyperbolic grid 
generation met hod in conjunction with an interactive grid generation system. The geometry and grid of the double- 
delta wing upper surface is shown in Fig. 17. The surface grid had 113 points in the chordwise direction and 129 
points in the circumferential direction w-ith an equal number of points on the upper and lower surfaces. This 
grid was clustered at the apex, the leading edges and the trailing edge by using hyperbolic tangent distributions. 
The leading edges were locally rounded to a diameter matching the wind tunnel model. A minimum of 13 points 
defined the curvature of the leading edges. The wing surface was used as the initial condition for the hyperbolic 
grid generation method to create the spherical C-0 type field grid, shown in Fig. 18. The field grid was generated 
by marching in the direction normal to the surface. Close to the surface 21-25 points were taken to sufficiently 
resolve the boundary layer and predict, flow separation. The y+ value of the first point, off of the surface was no 
larger than 1 for the Reynolds numbers simulated. The outer boundary of the grid was taken 2-2.5 chord lengths 
away from the surface in all directions and 73 points were used in the normal direction to extend the grid to the 
far field. To further resolve the convection of vortical flow, a grid of the same family of points, with twice the 
number of points in each direction, was embedded into the major grid above the wing using the Chimera or Overset 
methodology (see inset of Fig. 18). The connectivity between the two grids was prescribed by using the Domain 
Connectivity Function Three-Dimensions (DCF3D) developed at the NASA Ames Research Center. The embedded 
grid extended from near the apex to the downstream outflow boundary in the chordwise direction, from close to 
the symmetry plane to the leading edges in the spanwise direction, and from the boundary layers edge (selected 
as the 21st. point off the surface) to a maximum of 0.6 chord lengths above the wing in the normal direction. The 
dimensions of the embedded grid were 255 x 177 x 73 and the composite grid totaled slightly over 3.4 million points. 

To solve the flow over the double-delta wing, the compressible thin-layer Navier-Stokes code OVERFLOW 
developed at the NASA Ames Research Center was used. The code has an option to use the diagonal form of the 
implicit. Beam and Warming type approximate factored scheme developed by Pulliam and Chaussee. OVERFLOW 
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was run at. .40,4 - 10 deg, at a Mach number of 0.2 and Re = 0.5 M with completely laminar flow. Each run took 
24 hours on a Cray C90 computer. Since the computational grid planes did not. line up with constant chord sections, 
the CFD solution was interpolated onto the spacial coordinates used in the experimental flowfield surveys. A code 
was developed, which used a stencil walk, an iterative search algorithm, and tri-linear interpolation to efficiently 
perform the interpolation. 

Pseudo-time stepping using time steps scaled to grid metrics and a CFL number maximum limit was used 
to accelerate the convergence of the simulation. For each of these cases, the L2 norm of the residuals dropped 
2 orders of magnitude by 1000 iterations, at which point the major vortical structures formed. The residuals 
remained constant for the additional 2000 iterations run, indicating that, there were small regions in the flowfield 
that remained unsteady. By the 1500th iteration, the secondary separation lines, caused by the strake and wing 
vortices had formed and were refined slightly until the final 3000th iteration. Inspection of the flowfield revealed 
that the portion of the flowfield outboard of the wing vortex secondary separation line changed slightly up to 
the last, iteration. These changes were caused by spurious accumulation and dissipation of regions of vorticity 
embedded within the wing leading edge feeding shear layer. As these regions of vorticity appeared above the wing 
surface, they perturbed the local boundary layer separation and surface pressure outboard of the wing vortex sec- 
ondary separation line. The strake vortex and wing vortex secondary separation lines, and local surface pressure 
inboard of the wing vortex secondary separation line, however, remained relatively constant with iteration. No 
such perceptable perturbations existed over the strake portion of the wing. Therefore, the solutions are considered 
converged and valid in all portions of the flowfield except in the region outboard of the wing vortex secondary sep- 
aration line. Time-accurate unsteady simulations would be required to capture the time evolution of these vort ical 
structures and their impact on the local flowfield as performed over a single delta wing by Gordnier and Visbal [23]. 


5 Comparison between Solution and Experiment 

Flowfield Data 

The vortex flow over the wing was surveyed at AO A - 10 deg and Re - 0.5 M in order to provide data for the 
validation of the Navier-Stokes solution at these conditions with fully laminar flow. The upper plots of Figs. 19 to 
21 show the velocity vectors, the contours of constant chord wise velocity component u/U^, and the total-pressure 
coefficient Cp t measured at 50, 62.5 and 75% chord. Cp t has been defined here as the difference between the local 
total pressure and the free stream static pressure, non-dimensionalized by the freestream dynamic pressure. The 
computed data is presented in the lower plots. The chordwise stations are located well upstream of the bursting 
point of the strake and wing vortices at this angle of attack. The laserlight-sheet flow visualization indicated that, 
the bursting point can be pushed upstream by the 5HP when the angle between the probe and vortex axes is large. 
This undesired effect was avoided by carrying out the surveys at 50, 62.5 and 75%> chord with the probe yawed at 
angles of 7, 20 and 15 deg, respectively. Since the rotation about the pitch axis was not computer controlled, the 
probe pitch angle was kept at zero deg during the surveys. At 75% chord, this unfortunately resulted in a void in 
the data measured beneath the wing vortex. In this region probe flow angles were encountered that were outside 
the calibration range of the 5HP. Due to the symmetry of the flowfield, the data is presented for the right half of 
the double-delta wing only. 

The vector plots illustrate the development of the vortex flow which exists downstream of the kink. For clarity, 
these plots do not show all the measured and computed data points. The development of the strake and wing 
vortices from 50 to 75% chord is reasonably well predicted by the solution. On the strake. the primary vortex 
core is predicted to be slightly outboard of the measured location, while on the wing the strake and wing vortex 
cores are predicted to be slightly inboard of the measured locations. The reason for this will be discussed in the 
subsequent section. 

The u/Uzc and Cp t measured in the strake vortex changes hardly from 50 to 62.5% chord. From 62. 5 to 75% 
chord, the magnitude of this data can be seen to decrease slightly. This may be due to the fact that at some 
downstream distance from the kink the strake vortex is no longer fed with vorticity from the strake leading edge. 
The wing vortex, on the other hand, increases rapidly in strength with downstream distance from the kink and at 
62.5% chord in this vortex much larger u/U cc and lower Cp t values are recorded than in the strake vortex. These 
tendencies are well predicted by the solution. The magnitude of the computed quantities is in general lower than 
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that measured in the experiment. 

Surface Flow 

Particle traces, resembling surface oil flow visualization, for the laminar flow solution at ,40.4 = 10 deg and 
Re = 0.5 M are compared with the experiment in Fig. 22. On the strake, the secondary separation line and the 
flow that is entrained is well predicted. A tertiary separation line is also predicted whose position agrees wbth 
the experiment. Flow entrainment towards the tertiary separation line in the outboard direction is weak due to a 
lack of grid resolution. At the kink, the surface flow resembles the experiment, in that the strake vortex tertiary 
separation line vanishes, and the flow is entrained outboard towards the wing vortex secondary separation line. A 
slight outboard bending of the strake vortex secondary separation line near the kink, caused by the influence of 
the formation of the wing vortex, is also predicted. 

Over the wing, the strake vortex secondary separation line continues to exist. The strake vortex secondary 
separation line is predicted to be a few percent semispan inboard of the experimental position. Downstream of the 
85% chordwise location, the predicted secondary separation line bends inboard to become parallel to the centerline 
at the trailing edge. The wing vortex secondary line is also predicted inboard of the experimental position. The 
outboard bending of the wing vortex secondary separation line noted in the experiment was caused by transition 
from laminar to turbulent flow. Since the flow simulated was purely laminar, this outboard bending could not 
be predicted. The position of the simulated wing vortex secondary separation line was therefore inboard of the 
experimental line, which forced the strake vortex and secondary separation line inboard as well. Another key 
feature of the surface flow captured is the apparent focus near the cropped portion of the wing. In this location, 
the solution is converged and flow issues from a point in all directions, but predominantly upstream, as observed 
in the experiment. Upstream and outboard of this focus, the surface flow prediction is unreliable due to local lack 
of convergence, as discussed in the section on Computations. 


Surface Pressure 

Fig. 23 compares the measured and computed Cp at 25, 75 and 90% chord for AO A — 10 deg and Re ~ 0.5 M. 
The shape of the Cp distributions induced by the strake primary and secondary vortices is in general well predicted 
by the solution. At 25% chord, the strake vortex suction peak is predicted to be far outboard of the measured 
location. It should be reminded here, that the effective AO A may be on the order of 12 deg in the experiment. A 
reduction of AO A from 12 to 10 deg, decreases the strake vortex suction and moves its peak outboard, but not so 
much that it matches the computed peak. At 75 and 90% chord, both the strake and the wing vortex are predicted 
to be slightly inboard of the location found in the experiment. A reduction of the experimental AO A from 12 to 10 
deg increases the Cp induced by the strake and wing vortices by about 0.2. The differences between the predicted 
and experimental Cp distributions hence become smaller when this is taken into account. 

Fig. 25 shows coloured pictures of the computed Cp distribution on the entire upper surface of the double-delta 
wing and the Cp distribution that was determined from PSP test results. The PSP results shown were obtained 
at Re = 1.5 M, because the resolution at this Re was better than that at lower Re. Since the computations were 
performed at Re = 0.5 M, a direct comparison with the PSP results is not justified. The PSP picture shows white 
spots on the central part of the strake, suggesting a positive Cp in this region. This is incorrect and due to glare 
from the lamps that were used to excite the PSP. Both pictures clearly show the negative-Cp regions induced by 
the strake and wing vortices. The wing vortex can be seen to induce a lower Cp than the strake vortex. In the 
PSP experiment, a Cp as low as -3.1 was induced by the wing vortex just downstream of the kink. This agrees well 
with data measured on an identical double-delta wing at TUD. where under similar test conditions a minimum Cp 
of -2.9 was recorded at 60% chord [22]. The computed wing vortex pressure field can be seen to suddenly expand 
at about the 85% chordwise position. This correlates well with the vortex core bursting location observed in the 
laserlight sheet flow visualization tests(Fig. 8). 


Turbulent Flow Computations 

Currently, computations are performed with OVERFLOW running at Re = 1.0 M with fully turbulent flow. To 
simulate completely turbulent flows, the Baldwin-Barth one equation turbulence model is used. The first results 
indicate that the cores of the strake and wing vortices are predicted to be outboard when compared to the laminar 
case. This is accompanied with an outboard shift of the secondary separation lines of both vortices. At 25% chord, 
the surface Cp is stronger negative for the turbulent case, while the Cp induced by the secondary vortex is less 


negative. This is in accordance with earlier experimental results [21]. At 75% chord, the Cp induced by the st.rake 
vortex is nearly identical for either the laminar or turbulent case (Fig. 24). The wing vortex suction, however, is 
stronger and more outboard compared to the laminar flow case. 


6 Conclusions and Recommendations 

A low-speed wind-tunnel study of the flow about a 76/40-deg double-delta wing is described that was carried out 
at angles of attack ranging from -10 to 25 deg and at Reynolds numbers from 0.5 to 1.5 M. Flow visualization tests 
have provided an insight into the effect of the angle of attack and Reynolds number on the vortex-dominated flow 
on and off the surface of the double-delta wing. Surface pressure recordings and PSP tests have provided data on 
the pressures induced by the vortices. Flowfield data were obtained at an angle of attack of 10 deg by surveying 
the flow over the model using a thin 5-hole probe. The data for this angle of attack has been used to validate the 
solution of the thin-layer Navier-Stokes equations for fully laminar flow. The solution is shown to well predict the 
behavior of the flow on and off the upper surface of the double-delta wing. 

For a more complete validation of the numerical methods, the experimental tunnel wall and model support 
interference effects should be quantified. The numerical methods, on the other hand, should be extended in order 
to be capable to calculate these effects. 

For further validation of the numerical solutions for the flow over the present double-delta wing, additional 
experimental data are required, e.g., surveys of the vortex flowfield at angles of attack other than 10 deg and 
balance measurements. 
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Figure 1: BART wind tunnel 



Figure 2: Model geometry (dimensions in mm) 



Figure 3: Model in wind tunnel 
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Figure 4: Vortex flow illuminated by laserlight sheet; AO A = 10 deg 
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Figure 9(a): Upper surface flow pattern; AO A = -10 deg 
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Figure 9(c): Continued; AO A — 10 deg 
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Figure 9(d): Concluded; AO A = 22.5 deg 
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irface flow pattern; .40,4 =10 and 22.5 deg 
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Figure 13(b): Continued 
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Figure 13(c): Concluded 
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Figure 15: Effect of angle of attack on upper surface Cp at 90% chord; Re = 1.0 M 
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Figure 16: Reynolds number effect, on upper surface Cp at 75% chord: .40/1 = 10 and 22.5 deg 
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Figure 17: Surface grid 



Figure 18: Field grid 
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Figure 19: Measured and computed velocity vectors at AO A — 10 deg 
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Figure 20: Measured and computed u/U 0 0 at ,40.4 = 10 deg 
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Figure 21: Measured and computed Cp t at AO A — 10 deg 
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Figure 24: Comparison of laminar and turbulent flow solution for upper surface Cp at AO A = 10 deg 
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Figure 25: Upper surface Cp at AO A = 10 deg; computation and PSP results 
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